Emergence of domains and nonlinear transport in the zero-resistance state 
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We study transport in the domain state, the so-called zero-resistance state, that emerges in a 
two-dimensional electron system in which the combined action of microwave radiation and magnetic 
field produces a negative absolute conductivity. We show that the voltage-biased system has a rich 
phase diagram in the system size and voltage plane, with second- and first-order transitions between 
the domain and homogeneous states for small and large voltages, respectively. We find the residual 
negative dissipative resistance in the stable domain state. 

PACS numbers: 73.50.Fq, 64.60.an, 73.50.-h, 64.60.-i 
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Introduction. — The zero-resistance state (ZRS) |l|-|4( is 
perhaps the most spectacular manifestation of the newly 
discovered nonequilibrium effects in ultrahigh mobility 
two-dimensional (2D) electron systems in high Landau 
levels, driven by ac (in the microwave range) or dc elec- 
tric fields [f| . ZRS is attributed @ to the instability of a 
homogeneous state with the negative absolute dissipative 
conductivity a < and the associated nonequilibrium 
phase transition into a static domain state with zero net 
resistance [7Hl"2j]. The domain picture is supported by a 
number of experiments 13-2(J. Similar electrical insta- 



bilities have also been known to appear in other contexts 
[E Hl| i most prominently, in the Gunn diode [22| (where 
most of the effects are due to the negative differential 
conductivity and the emergence of moving domains) and 
in illuminated ruby crystals 23, l24| (where the strongly 
anisotropic nature of charge transport reduces the prob- 
lem to one dimension, with the instability controlled by 
the differential conductivity). What also makes the ZRS 
problem special is the emergence of a < in the presence 
of a magnetic field B that produces a strong Hall compo- 
nent of current. While the microscopic mechanisms that 
lead to (J < in nonequilibrium 2D electron gases at 
B 7^ are by now fairly well established the physics 
of the resulting domain state remains poorly understood. 

Most works so far have studied the bulk proper- 

ties of the domain state, i.e., the limit d/L — > in which 
the width d of the domain wall (DW) is vanishingly small 
compared to the system size L. In fact, however, the 
physics related to the DW structure is crucially impor- 
tant near the phase transition, because the DW width 
diverges at the critical point and serves as a critical pa- 
rameter of the transition. Moreover, the position of the 
phase transition in finite-size systems (nonzero d/L) be- 
comes L dependent and shifts towards negative a [251 ] . 

In this paper, we develop an analytical model of the 
domain state for arbitrary d/L and study the nonlinear 



response of the domain state to external voltage. The 
analytical solution enables us to construct the phase dia- 
gram of the biased finite-size system, which incorporates 
not only continuous but also discontinuous transitions 
between the homogeneous and domain states, and to cal- 
culate the negative conductance in the domain state. 
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FIG. 1. (a) Geometry of the model: a 2D stripe between two 
metallic contacts along the long sides. The double dashed 
line shows a possible position of the domain wall, the arrows 
indicate the direction of the electric field E and the Hall cur- 
rent ]h in the voltage-biased domain state, (b) The negative 
nonlinear dissipative conductivity o~(E) and the diffusion co- 
efficient D independent of E. 

Model. — Consider 2D electrons occupying a stripe (\x\ < 
L/2, z = 0), infinite in the y direction (Fig.QJi), in be- 
tween two plane metallic contacts at x = ±L/2 that are 
perpendicular to the stripe. Under the illumination by 
microwaves at B ^ 0, the linear response dc dissipative 
conductivity a(E — > 0), where E is the dc driving field, 
becomes negative in one or more intervals of B for the mi- 
crowave power P above the critical threshold P c [H[ . The 
nonlinear absolute conductivity <j(E) remains negative 
in a finite range of E, crossing zero at the critical field 
E c (Fig. [lb) [26|. By contrast, the diffusion coefficient 
D > is nearly unaffected by the radiation and obeys 
the Einstein relation <7dark = g 2 \D valid at equilibrium, 
where Odark and \ = m/it (m is the electron mass) are 



2 



the dark conductivity and the dark compressibility. 

We first explore the part of the phase diagram in which 
the system at P =^ remains homogeneous along the 
y axis, so that the surface electron density n e {x) and 
the x and y components of the electric current j x , y {x) 
depend on x only. The domain state is a stable solution of 
the Poisson and continuity equations. The former relates 
n e (x) to the normal component £ z {x,z) of the electric 
field at z = ±0. Namely, 2-Ke[n e (x) — no] = e£ z (x, +0) 
with e < 0, where no is the density of background positive 
charges and e is the dielectric constant of the medium. 
The continuity equation reduces in the static limit to 
d x j x (x) — 0. We assume that the in-plane component 
of the electric field E(x) = S x (x, +0) varies on a spatial 
scale that is larger than the microscopic scales [13] of the 
problem and introduce the local conductivity a[E(x)]. 
Importantly, at nonzero P, the Einstein relation does not 
hold and the sum of the drift and diffusion terms in the 
dissipative current j x =a(E)E—eDd x n e is not ex pres sible 
as the gradient of the electrochemical potential [28j. 

At z ^ 0, the functions £ x {x, z) and £ z (x,z) are har- 
monic conjugates. This, together with the continuity of 
£ x (x, z) at z = 0, allows us to represent the Poisson 
equation as eE{x) = 2i:e'H{n e (x)}, where the Hilbert 
transform H{f(x)} = 7r _1 p.v.J dx'(x — x')~ 1 f(x') obeys 
Ti 2 = —1. Applying the Hilbert transform to the Pois- 
son equation in this form and substituting the result in 
the diffusion term in j x , we find that E(x) satisfies the 
nonlinear integral equation 



E{x)a[E(x)} + (eD/2n) d x H{E{x)} = j x 



(1) 



Below, we solve Eq. ((T|) with the boundary conditions 
n e (x = ±L/2) = no and for a particular choice of 



a(E) = a{0)(E c /7rE)sin(irE/E c ) 



(2) 



with cr(0) < (plotted in Fig.Q})) |29j |. 
Domain solution. — By introducing the dimensionless 
density p(x) — ir 2 e{n e (x) — n ]/eE c , electric field 6{x) = 
itE(x)/2E c , and current j = Trj x /a(0)E c , Eqs. (1) and 
(2) are rewritten in the compact form as 

sm29 + 2Xd x p = j , (3) 

where p(x) and 9(x) satisfy p(x) = —H{9(x)} and 

A = eL>/27r|cr(0)| . (4) 

Being the only spatial scale in Eq. j3]) at j = 0, A is 
identified as the nonequilibrium screening length which 
gives the characteristic width of the DW [271 ] . 

With the boundary conditions p(±L/2) = 0, the do- 
main solution to Eq. ([3]) in terms of the complex function 
*dom = #dom + ipdam reads 



i In 



cosh(£ — iw) arcsinj 
sinh£ 2 



(5) 



where 2£ = inx/L + iw + (3, w = arctan(jZ), (3 = 
arcoth\A 2 — j 2 l 2 , and / = L/ttX. The system sponta- 
neously chooses between two degenerate domain states 
related to each other by $d om (i) <-» fj om (- x). For def- 
initeness, we analyze the solution with 9dam(L/2) > 0. 
Homogeneous state: Linear stability. — For given j in the 
homogeneous case of p — 0, the electric field 0hom satisfies 
sin 2#hom = j [Eq. ©] with the solution 



? hom = (V 2 ) arcsin j 



?Ln = */2 - t , (6) 



where the signs ^ correspond to the regions with the 
negative (<) and positive (>) differential conductivity 



a d (E) = d E [Ea(E)} = ct(0)cos26» . 



(7) 



Linear stability analysis [30| in the vicinity of this solu- 
tion (with the electric field along the x axis) shows that it 
is stable against small time- and space-dependent charge 
fluctuations proportional to exp(iq x x + iq y y) if 



a d (E)q 2 x +a(E)q 2 >-(eD/27r)(q 2 x 



^) 3/2 



(8) 



for all possible q x and q y . In an infinite 2D system, Eq. ([8|) 
reduces to the usual stability conditions ad > for lon- 
gitudinal (q y = 0, q x 0) and a > for transverse 
(q x = 0, q y —> 0) fluctuations. In the stripe geometry, 
q x takes discrete values im/L with |n| = 1,2,... and the 
diffusion term on the right-hand side of Eq. ([8]) becomes 
relevant. Equation ([8]) then yields 



a d {E) > -eD/2L 



(9) 



for the condition of the longitudinal stability {q y = 0, 
\q x \ = tt/X). For E — > 0, Eq. © gives the threshold value 
of / = 1 [25j for the breakup of the homogeneous state 
into domains in the unbiased discussed below. 

Unbiased domain state. — For j = 0, Eq. ([SJ reduces to 



#dom = arctan [(I 2 
Pdom = artanh [(1 



- l)^ 2 sin(Trx/L)} , 
P 2 ) 1/2 cos(7ra;/L)] 



(10) 
(11) 



(Fig.[2ji,b). In the limit I > 1, Eq. (10]) simplifies to 
$dom — arctan(x/A), which means two domains with 
E(±L/2) ~ ±E C separated by the DW of width A. The 
DW is charged with p{x) ~ \n(L/\x\) for A < \x\ < L. 
This gives the oppositely directed Hall currents j y (x) = 
-en e (x)cE(x)/B on the sides of the DW (Fig.[T]a,). With 
lowering I, both \8{x)\ and p(x) decrease, and vanish to 
zero at I = 1. As follows from Eq. (j9]), the homogeneous 
state with cr(0) < in the unbiased stripe is stable for 
I < 1, i.e., there is a continuous transition between the 
homogeneous and domain states. Near the transition, for 
< I — 1 1, ^dom — i\/2(l — 1) cxp(— inx/L) vanishes 
with the critical exponent 1/2. 

Biased domain state. — For j ^ 0, Eq. ([5]) tells us that the 
DW shifts by Lw/ir from x = while the characteristic 
width of the DW d = Lj3/w grows with increasing |j| 
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FIG. 2. Spatial distribution of (a) the electric field E(x) 
(in units of E c ) and (b) the charge density p(x) (in units 
of iv 2 e/eE c ) in the domain state of the unbiased stripe (cur- 
rent j x = 0) for L/tv\ — 30, 10, 2, 1.1, 1.01 (L decreases in 
the direction of arrow). As the current is increased, the do- 
main wall is shifted and broadened. The field E(x) is shown 
for (c) L/tt\ — 3 and the current [in units of a(Q)E c /Tr] 
j = 0, 0.25, 0.5, 0.75, 0.9, (8/9) 1/2 and for (d) L/ttA = 30 
and ] = 0, 0.03, 0.1, 0.5, 0.9, (899/900) 1/2 (j grows in the 
direction of arrow). 



(Fig.[2t,d) and diverges as (j C 2 
point \j\= j c2 , where 



\j\) 1 / 2 at the critical 



jc2 = (1 - I' 2 ) 



aa/a 



(12) 



According to Eqs. (0 and ([9]), for \j\ > j C 2 the homo- 
geneous state 0hom is stable against longitudinal fluctua- 
tions. The line (fT2"j) of the second-order phase transitions 
includes, at its endpoint at j c % — 0, the transition at 
I = 1 in the unbiased stripe discussed above. 

Averaging the field ^dom [Eq. ©] for given j over the 

stripe cross-section, 6>dom = {l/L)f_ L ^ 2 dx0dom{x), one 
finds the bias voltage V = 2E c L9d onl /Tr and the current- 
voltage characteristic (CVC) of the domain state with 



arctan(j7) — (1/2) arcsinj 



(13) 



Note that the current j x flows against the applied field, 
Jx V < 0. For V -»• 0, Eq. (USD gives 9 dom = (I - 1/2)], 
or, restoring units, the linear dissipative conductance of 
the stripe j x /V = {o~)/L in the domain state, where 



<r(0) 



2L/tt\- 1 



L > 7rA 



(14) 



see inset in Fig.[3J It is worth noting that, as L increases, 
(a) in Eq. (Tl4)) behaves as L^ 1 in sharp contrast to the 



exponential behavior of (a) cx — exp(— L/A3D) [8( for a 
three-dimensional medium with the negative conductiv- 
ity 03D i where the relation between the electric field and 
charge density is local. Here A3D = (eD/47r|(T3D|) 1 ^ 2 is 
the analogue of A [Eq. (Q}]. Transport across the DW is 
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FIG. 3. Phase diagram of the voltage-biased stripe in the V— 
L plane (V and L measured in units of E C L and 7rA, respec- 
tively). The phase boundary (solid line), induced by the lon- 
gitudinal instability, separates the homogeneous (above) and 
domain (below) states. In the unbiased stripe (V — » 0), there 
is a second-order phase transition at L/ttX = 1. The continu- 
ous transition line terminates with increasing V and L at the 
tricritical point (filled circle) V = E c L/4, L/tv\ — \f2. At 
larger V and L, the transition is first order. The shaded area 
is the region of hysteresis. The thin line, with the endpoint 
marked by the triangle, denotes the linear stability thresh- 
old for the homogeneous state against transverse fluctuations 
(the cross marks the point at which the line intersects the 
discontinuous transition line obtained in the model of frozen 
transverse fluctuations). Inset: the effective linear-response 
conductivity (a) as a function of L [Eq. (|14|l ]. For compari- 
son, the dashed line shows the behavior for a 3D model 



thus strongly enhanced by the nonlocal character of 2D 
electrostatics. 

The CVC (fT3"]) for several values of / is shown by dashed 
lines in Fig.UJ For I > y/2, the current is seen to become, 
as the voltage is increased, a double-valued function of 
V. That is, in fact, the continuous transition line (fT2")) 
in the V-l plane terminates at I = \[2 and becomes first 
order for larger I (Fig. 13]), as we discuss next. 

Lyapunov functional. — The linear stability analysis, 
which leads to Eq. ^ and the continuous transitions 
in the interval 1 < I < y/2, does not capture the emer- 
gence of the discontinuous transitions for / > y/2, i.e., [by 
substituting I — \/2 in Eqs. (fl~2)) and (fT3))] for large volt- 
ages V > E C L/A. The stability analysis of the domain 
solution ([5]) that we perform below to describe the large- 
voltage regime is based on the Lyapunov functional (LF) 
approach to the ZRS problem [llj. The advantage of 
the LF method is that it is capable of discriminating the 
stable (global minimum of the LF) and metastable (local 
minimum) states that can be distant in phase space. The 
LF ${_E(x)} = — G + K is given by the difference of the 



gain G ■■ 
tion K 



-Jdx J^ {x) dE'E'a(E') and the DW contribu- 
(D/2) Jdx E(x)CE(x), where the capacitance 
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FIG. 4. Current-voltage characteristic (thick lines) of the 
stripe for different L/ir\. For L/ttX < 1 (homogeneous state), 
the dependence of the dissipative current jx{V) is given by 
the lowest curve j = smfjrV/EcL). For 1 < L/tyX < \/2, 
the homogeneous state breaks up into domains for V < V C 2, 
in a continuous fashion, and the j x —V curve has a kink, as 
shown for L/ttX — 1.2. For L/nX > y/2, the j x — V curves 
for the domain state (dashed lines) are double-valued and 
the transition becomes discontinuous, as demonstrated by the 
jumps of jx(V) at V = V c i (the critical voltage is marked for 
L/irX = 10) in the curves for L/ttX = 2.5, 10, 100, 1000. The 
discontinuous transition line (dash-dotted) terminates at the 
tricritical point (filled circle) at V = V C 3- The arrows on the 
thin vertical lines denote hysteresis (shown for L/nX — 10) in 
the interval V m i a < V < V max . 



operator C = (e/2Tr)d x H. For the model ([2]), we have 

~6{x)-6(x') 



' dx , o n , s . f dxdx 1 



, (15) 



which gives $ho 



$dom = — I 1 



sin (9 nom in the uniform state and 



I-P )/2+l- 1 -l- 1 lnJf + l-^ (16) 



in the domain state §5§ [33 ]. 

First-order transitions. — Comparison of $hom and $dom 
leads to the phase boundary (thick line) in Fig. [3] and the 
CVC (thick lines) for different / in Fig. H For / < 1, the 
homogeneous state is stable ($hom < ^dom) for arbitrary 
V. For a given I in the interval 1 < I < v2 f there is 
a continuous voltage-driven transition ($hom = ^dom) 
between the homogeneous and domain states at V = 
V C 2 = (E c L/ir) arccos/ -1 . For I > y/2, there emerges 
the interval Vmi n < V < Vmax (whose endpoints are 
shown in Fig. 0] for I — 10) in which the function j(V) 
for the domain state is double-valued (dashed lines in 
Fig HJ- Note that the expression for Vmin coincides with 
that for V C 2, so that V m - m saturates as I is increased at 
E c L/2. On the lower branch (with dV/dj < 0), $dom 
is larger than on the upper (and the lower branch is 
unstable against linear fluctuations [30|). For the up- 
per branch, the phase boundary equation <I>hc 



(whose solution is shown in Fig. [4] as a dash-dotted line) 
yields, for I > y/2, the first-order transition at V = V c i, 
where V c \ tends to E C L in the "bulk" limit I — > oo as 
V cl /E C L ~ 1 - (l/7r)(21nZ/Z) 1/2 [33]. The current dis- 
continuity [which vanishes for j as (21n//Z) 1 / 2 for I ^> 1], 
between the upper-branch domain state and the homo- 
geneous state, is illustrated by the vertical thick lines 
in Fig. |U The "tricritical" point, separating the first- 
and second-order transitions, at / = v2 and V — V C 3 = 
E c L/4 is marked by the filled circle in Figs. [3] and 01 

The domain state for V c i < V < V msx and the ho- 
mogeneous state for Vmi n < V < V c i are metastable, i.e., 
they can be probed if the voltage sweep rate is larger than 
their characteristic decay rates (the dynamical properties 
of the system are beyond the scope here: in the above, 
we assumed that the system resides in the stable state, 
in which $ is minimized globally). In the nonadiabatic 
limit, the system exhibits hysteresis [3l|, as marked in 
Fig. |4] by the arrows on the thin vertical lines for the 
case of I = 10. The hysteresis range V" m i n < V < V m ax 
for arbitrary I > v2 is shown as a shaded area in Fig. [3] 
Transverse instability. — Before concluding, we briefly 
comment on the stability of the above picture against 
transverse fluctuations. Allowing for nonzero q y in 
Eq. ((SJ , they can be shown [30] to be irrelevant on the 
phase boundary (thick line) in Fig. [3] for a sufficiently 
narrow stripe, namely for I < l c ± ~ 1.76 (or, equiva- 
lently, V < V c ± ~ 0. 31 E C L). The threshold is marked in 
Fig. [3] by the triangle. The range of I < l c ± includes the 
zero-bias critical point at I = 1 and the tricritical point 
at I = v2- For / > l c ±, the linear transverse stability is 
maintained above the thin line in Fig. G2 which runs well 
above the longitudinal stability threshold Vmin for the 
homogeneous state ([9]) and very closely to the discon- 
tinuous transition line V c i obtained for frozen transverse 
fluctuations. At I = l c+ ~ 7 (marked by the cross in 
Fig-EJ), the two lines intersect, so that for I > l c+ the ho- 
mogeneous state is unstable above the phase boundary 
at V c i. In the narrow region between the lines for I > l c+ , 
the global minimum of the LF should thus be given by 
a 2D domain state with broken translational invariancc 
along the stripe. The nature of this state, as well as the 
position of the boundary V c \ for the global stability of 
the domain state 34|, requires additional study. 
Summary. — We have studied transport in the voltage- 
biased stripe with a negative absolute conductivity and 
obtained the phase diagram which shows phase transi- 
tions between the domain and homogeneous states. The 
transitions are second order for small and first order for 
large voltages (Fig. [3]). We have calculated the current- 
voltage characteristic of the domain state (Fig. [4]) and 
found the negative dissipative conductance. Our predic- 
tions can be verified by measuring the current-voltage 
characteristic in sufficiently small samples in the vicinity 
of the ZRS transition. 

Acknowledgements. — We thank S. Dorozhkin, 



5 



Y. Galperin, A. Kamenev, P. Ostrovsky, I. Pro- 
topopov, J. Smet, R. Suris, and M. Zudov for interesting 
discussions. The work was supported by the DFG- 
RFBR (I.A.D. and A.D.M.) and by the University of 
Iowa (M.K.). 



[1] R.G. Mani, J.H. Smet, K. von Klitzing, V. Narayana- 

murti, W.B. Johnson, and V. Umansky, Nature (London) 

420, 646 (2002). 
[2] M.A. Zudov, R.R. Du, L.N. Pfeiffer, and K.W. West, 

Phys. Rev. Lett. 90, 046807 (2003). 
[3] C.L. Yang, M.A. Zudov, T.A. Knuuttila, R.R. Du, L.N. 

Pfeiffer, and K.W. West, Phys. Rev. Lett. 91, 096803 

(2003) . 

[4] D. Konstantinov and K. Kono, Phys. Rev. Lett. 105, 
226801 (2010). 

[5] LA. Dmitriev, A.D. Mirlin, D.G. Polyakov, and M.A. 

Zudov, Rev. Mod. Phys. 84, 1709 (2012). 
[6] A.V. Andreev, LL. Aleiner, and A.J. Millis, Phys. Rev. 

Lett. 91, 056803 (2003). 
[7] M.G Vavilov and LL. Aleiner, Phys. Rev. B 69, 035303 

(2004) . 

[8] A.F. Volkov and V.V. Pavlovskii, Phys. Rev. B 69, 
125305 (2004). 

[9] J. Alicea, L. Balents, M.P.A. Fisher, A. Paramekanti, and 

L. Radzihovsky, Phys. Rev. B 71, 235322 (2005). 
[10] LA. Dmitriev, M.G. Vavilov, LL. Aleiner, A.D. Mirlin, 

and D.G. Polyakov, Phys. Rev. B 71, 115316 (2005). 
[11] I. Finkler, B.I. Halperin, A. Auerbach, and Y. Yacoby, J. 

Stat. Phys. 125, 1093 (2006). 
[12] LA. Finkler and B.I. Halperin, Phys. Rev. B 79, 085315 

(2009). 

[13] R.L. Willett, L.N. Pfeiffer, and K.W. West, Phys. Rev. 

Lett. 93, 026804 (2004). 
[14] M.A. Zudov, R.R. Du, L.N. Pfeiffer, and K.W. West, 

Phys. Rev. Lett. 96, 236804 (2006). 
[15] W. Zhang, M.A. Zudov, L.N. Pfeiffer, and K.W. West, 

Phys. Rev. Lett. 98, 106804 (2007); A.T. Hatke, H.- 

S. Chiang, M.A. Zudov, L.N. Pfeiffer, and K.W. West, 

Phys. Rev. B 77, 201304(R) (2008). 
[16] S. Wiedmann, GM. Gusev, O.E. Raichev, A.K. Bakarov, 

and J.C. Portal, Phys. Rev. Lett. 105, 026804 (2010). 
[17] A. A. Bykov, JETP Lett. 91, 361 (2010). 
[18] S.I. Dorozhkin, L. Pfeiffer, K. West, K. von Klitzing, and 

J.H. Smet, Nature Phys. 7, 336 (2011). 
[19] D. Konstantinov, A. Chepelianskii, and K. Kono, J. 

Phys. Soc. Jpn. 81, 093601 (2012). 
[20] A.T. Hatke, M.A. Zudov, J.D. Watson, and M.J. Manfra, 

Phys. Rev. B 85, 121306(R) (2012). 
[21] V.L. Bonch-Bruevich, LP. Zvyagin, and A.G Mironov, 

Domain electrical instabilities in semiconductors (Con- 
sultants Bureau, New York, 1975). 
[22] J.B. Gunn, Solid State Commun. 1, 88 (1963). 



[23] P.F. Liao, A.M. Glass, and L.M. Humphrey, Phys. Rev. 

B, 22, 2276 (1980); S.A. Basun and A. A. Kaplyanskii 

and S.P. Feofilov, JETP Lett. 37, 586 (1983). 
[24] M.I. Dyakonov, JETP Lett. 39, 185 (1984); M.I. 

Dyakonov and A.S. Furman, Sov. Phys. JETP 60, 1191 

(1984). 

[25] S.I. Dorozhkin, LA. Dmitriev, A.D. Mirlin, Phys. Rev. B 
84, 125448 (2011). 

[26] For P - P c < P c , the field E c scales as |cr(0)| 1/2 [|. 

[27] Using the local conductivity is justified if the nonequi- 
librium screening length A [Eq. (Q] is larger than the 
microscopic scales (the cyclotron radius R c and the in- 
elastic scattering length k n [f|). Note that A diverges at 
<j(0) — > 0, thus serving as the critical parameter for the 
transition to ZRS in the bulk limit L/X — > oo. Deep in the 
domain state, for |cr(0)| ~ o"dark, the screening length is 
of the order of Adark = -Mcr(o)=cr d k = e/2me 2 . In typical 
experimental conditions, Adark <C kn,Rc and the width 
of the domain wall is determined by l- m and R c - 

[28] LA. Dmitriev, S.I. Dorozhkin, and A.D. Mirlin, 
Phys. Rev. B 80, 125418 (2009). 

[29] This model function allows for an exact analytical solu- 
tion to Eq. (1), while the physics the solution reveals is in- 
sensitive to the exact shape of the single-parameter func- 
tion cr(E). Note that the behavior of cr(E) for \E\ > E c 
plays no role in what follows. 

[30] LA. Dmitriev, M. Khodas, A.D. Mirlin, and D.G. 
Polyakov, unpublished. 

[31] Note that the obtained hysteresis bears close similarity to 
the noncquilibrium first-order phase transition proposed 
in Ref. [24| to describe the experiment [53] on illuminated 
ruby crystals. The key differences between our problem 
and that in Ref. [2J| are the importance of transverse 
fluctuations, which are absent in the model of Ref. [241 ]. 
and the nonlocal character of 2D electrostatics. We also 
found a richer behavior of the system by studying the case 
of smaller system sizes, L/irX ~ 1, where there emerge 
the continuous phase transitions and the tricritical point 
that separates the regimes of the first- and second-order 
phase transitions. 

[32] Note that the absolute minimum <& = — 1 of the LF (p~5|) 
is realized in the state \0(x)\ = n/2, i.e., in a system of 
size L, only at \V\ = E C L. 

[33] It is worth emphasizing that our consideration was per- 
formed in the mean-field framework. The impact of fluc- 
tuations on this problem remains to be studied. In par- 
ticular, one can apply the methods of Ref. [f| where it 
was found (in the limit V — > and L — > oo) that fluc- 
tuations may force the ZRS transition to be first order. 
Remarkably, we obtain the first-order transition driven 
by finite voltage already on the mean-field level. 

[34] We were able to check numerically that the linear stabil- 
ity threshold of the domain state ([5]) remains above V c i, 
i.e., the state ([5]) is linearly stable in the whole region 
V < Vci, at least for I < 40. 



